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The strength functions of the double giant dipole resonances (DGDR) in 16 O and 40 Ca arc 
calculated with the use of an extended version of the time-dependent Hartree-Fock theory known as 
the time-dependent density-matrix theory. The calculations are done in a self-consistent manner, in 
which the same Skyrme force as that used for a mean-field potential is used as an effective interaction 
for a two-body correlation function. It is found that the DGDR in le O has a large width due to 
the Landau damping, although the centroid energy of the strength distribution is close to twice 
the energy of the giant dipole resonance (GDR) calculated in RPA. The DGDR in 40 Ca is found 
more harmonic than that in 16 O: the strength function of the DGDR in 40 Ca is similar to what is 
predicted from the strength function of the GDR in RPA. 
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The double phonon states of giant resonances have become the subject of a number of recent experimental and 
theoretical investigations | l],^] . Microscopic calculations of the strength functions of double giant dipole resonances 
(DGDR) have also been done based on the shell model and quasiparticle-phonon models However, few 

microscopic studies have been reported, in which a single-particle basis and a residual interaction are treated in such a 
self-consistent manner as used in random-phase-approximation (RPA) calculations for giant resonances ||. We have 
recently proposed a self-consistent approach based on an extended version of the time-dependent Hartree-Fock 
theory (TDHF) known as the time-dependent density-matrix theory (TDDM) in which the same Skyrme force 
as that used for the calculation of a mean-field potential is used as a residual interaction for a two-body correlation 
function. We applied the model to the double giant quadrupole resonances (DGQR) in 16 and 40 Ca || and showed 
that the DGQR's in these nuclei have strong harmonic properties. The aim of this paper is to report the results of 
the application of the TDDM approach to the DGDR's in 16 and 40 Ca. 

The formulation of TDDM is based on the truncation of the hierarchy of reduced density matrices, in which genuine 
correlated parts in a three-body density matrix and higher reduced density matrices are neglected [|l0| . The TDDM 
equations thus determine the time evolution of a one-body density matrix p and a two-body correlation function Ci 
defined by C% — pi — ^[pp] 7 where ^[pp] is the antisymmetrized product of the one-body density matrices and is 
a two-body density matrix. In TDDM, further truncation is made by expanding p and C2 with a finite number of 
single-particle states {ip a } as 

p(ll',<) = ^n aa ,(i)VaM)V&(l',i), (1) 



C 2 (12l'2',t)= ]T C Q ^^(*)^(l,t)#(2,t)C'(l'^)^(2',t), (2) 

a/3a'f3' 

where the numbers denote space, spin and isospin coordinates. The time evolution of p and C2 is determined by the 
following three coupled equations |8y: 

in— ^ a (i,t) = MMM*(i,t), (3) 

at 

ihh aa , = y^J(aP\v\'y6}C lSa >/3 - C a f3 7 s(jS\v\a' (3)}, (4) 
0-yS 
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ihC a p a '0i = B a p a ipi + Patpot'P' + H a p a >pi, (5) 

where h(l,t) is the mean-field hamiltonian and v the residual interaction. The term B a /3 a ipi on the right-hand side 
of Eq.(5) represents the Born terms (the first-order terms of v). The terms P a /3a'f3' and H a p a :pi in Eq.(5) contain 
C a /3a'f3' and represent higher-order particle-particle (and hole- hole) and particle- hole type correlations, respectively. 
Thus full two-body correlations including those induced by the Pauli exclusion principle are taken into account in the 
equation of motion for C Qj a Q '^'. The explicit expressions for B a p a /pi, P a f3 a '/3' and Hafia 1 p 1 are given in Ref. Q. The 
small amplitude limit of TDDM was investigated jll| and it was shown that TDDM can be reduced to the second 
RPA [l2]| - |]l4| in such a limit. The number of two-body matrices treated in TDDM grows very rapidly with increasing 
mass number, restricting the application of TDDM to light nuclei for the present. To solve the coupled equations, we 
use the Skyrme interaction of the form jl^] 

v(r - r') = t (l + x P a )S 3 (r - r') + ^{fc'^r - r') + <5 3 (r - r')k 2 } 

+ t 2 k' • 5 3 (r - r')k + \u P {*-^-\ S 3 (r - r'), (6) 

where k = (V r — V T ')/2i acts to the right and k' = (V r ' — V r )/2i acts to the left. The factor 1/2 on the density 
dependent term contains the contribution of a rearrangement effect jl5| . We use the parameter set of the Skyrme 
III force (SKIII) [ fL6| . The spin-orbit force is neglected. We assume that the motion of the DGDR is generated by a 
two-body operator D 2 : 

|*(t = 0))= e ^ 2 |$o), (7) 

where D is a one-body dipole operator and |$o) the ground-state wave function. The initial conditions for solving the 
coupled equations Eqs.(3)-(5) are determined with the use of the above wave function. We evaluate the initial values 
of C a p a >p* 

C a p a .p.{t - 0) - (*(t = 0)\a+ a%a p a a \*(t = 0)), (8) 

assuming that |$o) is the Hartree-Fock (HF) ground-state wave function. At first order of k, the initial condition for 
C a f3a'(3' becomes 

C^ U pa = (^la+a+a^a^l*) 

= 2ik{{(i\D\p)(v\D\a} - (fj.\D\(T)(v\D\p)} (9) 
Cp ali u = (vPla+a+acrapl*) 

= -2ik{(p\D\n)(a\D\p) - ( p \D\v)(<t\D\vl)}, (10) 

where p and a refer to unoccupied single-particle states, and p and v refer to occupied ones. We choose the dipole 
operator in the above equations such as D — t z z. Other elements of the initial C a f3 a '/3' vanish at first order of k. 
Similarly, non-varnishing initial values of n aa i become 

n w = {^\a+a^) 

= 2ikJ2(»\ D W){"\ D \p) ( n ) 

n Pt . = (*\a+a p \V) 

= -2ikY,{p\D\y){y\D\p). (12) 

For the initial ip^s we use the HF single-particle wave functions. The strength function of the DGDR, defined by 

S 2 (E) = \{^n\D 2 \<P )\ 2 S(E - En), (13) 

n 

is given by the Fourier transform of the time-dependent two-body dipole moment 1^2 (£) as 

1 P°° Ft 
S 2 (E) = — D 2 (t)sm—dt, (14) 
irkn J n 
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where D2 is given by 

D 2 (t)^{^(t)\D 2 \^(t)) 

= Y J {a\D 2 \a')n a , a + £ (a\D\a')(fi\D\p){A[n a , a np. f )] + C a ,p a p}. (15) 

act' a{3a' 0' 

The terms without C a p a i pi in the above equation have negligible contribution to the Fourier transformation in Eq.(14). 
The k dependence of 82(E) thus obtained is negligible as long as k is sufficiently small. The energy- weighted sum 
rule (EWSR) for the DGDR is given as 

f°° 1 
J ES2(E)dE = -($o\[D 2 ,[H,D 2 }\<P ) 

-($o|^ 2 |*o) + 4(ti + t2)(<S>o\RD 2 \<S>o), (16) 
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where H is the total hamiltonian and m the nucleon mass. The second term on the right-hand side of Eq.(16) is due 
to the momentum dependence of the Skyrme force and R is the following two-body operator 

R= Yl * 3 (n-*i)- ( 17 ) 

The EWSR value is evaluated with the use of the HF wave function for |$o}- The second term on the right-hand side 
of Eq.(16) has a contribution of about 30% to the total EWSR value. 

To solve the coupled equations Eqs.(3)-(5), we use a minimum number of single-particle states: the Is, lp, 2s and Id 
single-particle orbits for 16 and the Is, lp, 2s, Id, 2p and If orbits for 40 Ca. To check the validity of such truncation 
of the single-particle space, we performed RPA calculations for the GDR strength functions in 16 O and 40 Ca using the 
time-dependent RPA equations || in the same truncated space. The obtained results were compared with those of the 
TDHF calculations which correspond to continuum RPA calculations Q . The fractions of the EWSR values depleted 
in the energy interval — 40MeV were turned out to be 90% in 16 O and 93% in 40 Ca, respectively. These values 
are sufficiently large and comparable with the TDHF values which are close to 100% [[17). However, the excitation 
energies of the GDR's in RPA were slightly larger than those in TDHF. To adjust the excitation energies of the 
GDR's to the TDHF values, we reduced the parameter x of the spin-dependent term of the Skyrme force Eq.(6). 
The obtained value of Xq was 0.3 instead of the original value of 0.45. We use this reduced value of Xq in the following 
calculations. The spin-dependent term of the Skyrme force has a negligible contribution to the mean-field potential in 
spin-isospin symmetric nuclei like 16 and 40 Ca considered here |l5j and, therefore, the single-particle wave functions 
are not affected by the reduction of the strength of the spin-dependent term. The integration in Eq.(14) is performed 
for a finite time interval of 1.5 ~ 2 x 10~ 21 s. As a result 82(E) has small fluctuations. To reduce the fluctuations in 
82(E), we multiply D 2 (t) by a damping factor e~ r */ 2 before performing the time integration. This corresponds to 
smoothing the strength function with a width T. We use T = lMeV. Other calculational details are explained in our 
previous publications |s|Jl8|] . 

The strength distribution of the DGQR in 16 calculated in TDDM is shown in Fig.[j] (thick solid line). The bump 
seen around E = 45MeV corresponds to the DGDR. The strength function Si (E) of the GDR obtained from the time- 
dependent RPA calculation is also shown in Fig|l] (dotted line) . The width of the GDR is small and nearly equal to the 
width due to the smoothing and the finite time integration as explained above. The thin vertical bar at E = 45.8MeV 
indicates the location of the DGDR strength predicted from the GDR shown in Fig.|l|. The fraction of the EWSR 
value of the DGDR depleted in the energy interval 10-60MeV is 82%. The centroid energy of the DGDR strength 
distribution between 35MeV and 55MeV is 44.6MeV. The energy difference AE = 45.8MeV-44.6MeV= 1.2MeV is 
small but slightly larger than the value AE k 0.8MeV obtained from the work done by de Souza Cruz and Weiss JlSf| 
using the generator coordinate method and the value AE < 0.7MeV obtained from the formula given by Bertsch and 
Feldmeire J2(|. Though the centroid energy of the DGDR is close to twice the GDR energy, the DGDR in 16 O has a 
large width due to the Landau damping. 

The strength distribution of the DGQR in 40 Ca calculated in TDDM is shown in FigJ| (solid line). The bump seen 
around E = 40MeV corresponds to the DGDR. The strength function Si(E) of the GDR in 40 Ca calculated in the 
time-dependent RPA is also shown in Fig.|| (dotted line). The GDR strength is split into two peaks in the case of 
40 Ca. A similar split is seen in the TDHF calculation for 40 Ca 0. The fraction of the EWSR value of the DGDR 
depleted in the energy interval 10-60MeV is 88%. In Fig.|| the strength function of the DGDR (thick solid line) is 
compared with what is expected from that of the GDR shown in Fig|2| The thin vertical bars in Fig.|| indicate the 
locations and relative strengths of the DGDR predicted from the GDR in RPA, assuming that the GDR consists 
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of the two discrete components. The centroid of the DGDR strength distribution in the energy range 30-50MeV is 
39.6MeV, while that of the GDR prediction is 39.0MeV. The energy difference AE = 0.6MeV may be larger than 



0.2MeV obtained from Refs. |19 20j] but is smaller than l.OMeV given by the shell model calculation for 40 Ca [|3|. 
Though the main peak in TDDM located at 40MeV has more strengths than that predicted from the RPA calculation, 
the DGDR seems to have weaker Landau damping in 40 Ca than in 16 O and the shape of the strength distribution is 
similar to what is expected from the GDR in RPA. This means that the dipole mode becomes fairly harmonic in the 
mass region of 40 Ca. A similar conclusion was obtained by de Souza Cruz and Weiss [[L9| comparing the DGDR in 
16 O with that in 40 Ca. The small anharmonicities of the DGDR in Ca have also been pointed out by Catara et al. 
pl[ using a boson expansion approach and the anharmonicities of the DGDR's in heavier nuclei have recently been 
studied by several groups |||||] J|f22|| . 

In summary, the strength functions of the DGDR's in 16 and 40 Ca were calculated in TDDM in a self-consistent 
manner, in which the same Skyrme force as that used for the calculation of the mean-field potential was used for 
the two-body correlation function. It was pointed out that the strength function of the DGDR is obtained from the 
Fourier transform of the time-dependent two-body dipole moment. It was found that in both nuclei the excitation 
energies of the DGDR's are very close to twice those of the GDR's. It was also found that the DGDR in 16 has a 
large width due to the Landau damping, indicating the anharmonicity of the dipole mode in 16 O. 
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FIG. 1. Strength function S 2 (E) of the DGDR in 16 calculated in TDDM (thick solid line). The dotted line denotes the 
strength function S\(E) of the GDR in RPA and the thin vertical bar at E — 45.8MeV indicates the location of the DGDR 
predicted from the GDR in RPA (in arbitrary units). 
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FIG. 3. Strength function 82(E) of the DGDR in 40 Ca calculated in TDDM (thick solid line) is compared with what is 
predicted from the strength function of the GDR shown in FigO (thin vertical bars in arbitrary units). 
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